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We compare two mechanisms for the dissociation of heavy mesons in an infinite 
quark-gluon plasma: dynamic Debye screening and multiple scattering. Using the 
uncertainty principle inspired by a Schrodinger-like equation, we find that the cri- 
terion Ob — — with Og// = oi{N(. + ■^) is parametrically true both for 
the dissociation of fast moving heavy mesons with a size as due to dynamic Debye 
screening as well as for mesons at rest in the medium. In contrast, we find that the 
criterion for the dissociation of heavy mesons due to uncorrelated multiple scatter- 
ing is parametrically ob ~ 1 — • Therefore, multiple scattering is a more 

efficient mechanism for the dissociation of heavy mesons in an infinite hot plasma. 
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I. INTRODUCTION 

The problem of dissociation of bound states in a hot QCD medium is of great importance 

in heavy ion colhsions as it provides evidence for the creation of the quark-gluon plasma in 

heavy ion collisions [1] . In order to get a better understanding of the properties of this state 

of matter it is necessary to establish criteria under which bound states are not allowed to 

exist or are broken apart. We focus on the study of two particular mechanisms which we 

believe to be the main causes of dissociation in the plasma: Debye screening and multiple 

scattering with constituents of the plasma, 
fl 

In Ref. [2|, the authors discussed J/il) suppression due to Debye screening by the quark- 
gluon plasma and the importance of this signature to diagnose quark-gluon plasma formation 
in heavy ion collisions. They found that the criterion for the dissociation of a J/ip at rest is 
rj/^ = 1.61r£) with r^i = the Debye screening length. In contrast, for a heavy quark- 
antiquark pair moving at velocity v in an infinite strongly coupled A/" = 4 SYM plasma, 
;he AdS/CFT calculation shows that the screening length goes like rc){v, T) ~ r£)(0, T)/ ^ 
3|. This suggests we may expect a criterion for the dissociation of heavy quarkonia in such 



br the dissociation due 



a strongly coupled plasma of the form Tqq ~ rj:,{^,T) j As 

to multiple scattering with the constituents of the plasma, Ref. ^| addresses the problem 
of heavy meson suppression in a finite dense QCD medium and predicts suppression of 
S-mesons comparable to that of D-mesons at transverse momenta as low as ~ 10 GeV. 

The calculation of the screening effect for mesons at rest is not enough to understand the 
suppression observed in the data from heavy ion collisions. It is also important to establish 
when a fast moving meson is broken apart due to the presence of the plasma. The latter case 
is the one explicitly addressed in this paper. For the screening effect, the full calculation is 
too complicated to be performed analytically and there is not an obvious way of getting a 
simple picture where we can get a good estimate based on the uncertainty principle as done 
in 21]. By considering the Dirac equation in light-cone coordinates and in light-cone gauge 
we propose a Schrodinger-like equation where we can rely on the uncertainty principle to get 
a sensible estimate for the criterion for dissociation of bound states. This approach requires 
the calculation of the effective field produced by a fast moving charge with respect to the 
plasma, which is done in the appropriate region of momenta where we keep terms only up to 
first order in ^ ~ 7"^. Even though in covariant gauge the effective field in the rest frame 



3 

of the charge moving relative to the plasma with f ~ 1 is highly anisotropic!^, we find that 
after a gauge transformation and going back to the rest frame of the plasma, the anisotropy 
is suppressed as inverse powers of 7 compared to that due to the Lorentz contraction. We 
find that the criterion — ^ is also parametrically true for fast moving heavy mesons. 

For multiple scattering we get a criterion for dissociation in terms of the saturation 
momentum of the system Q^, namely ~ l/Qs- For this purpose we first identify the 
typical time between interactions inside the meson r^. Then the criterion for dissociation 
is given by the statement that if the quarks inside the meson pick up enough transverse 
momentum during that period then the meson breaks up. We take as the natural scale 
for transverse momentum inside of the heavy meson l/a^- The transverse momentum 
broadening is given by Qs where the role of the lenght of the plasma is played by tb (assuming 
an infinite plasma). In the case of uncorrelated multiple scattering, the criterion Qg ~ l/a^ 
gives 

[7«e//ln^]5T' ^ 

with a^ff = a{Nc + ^). It is also parametrically true for the dissociation of heavy mesons 
almost at rest with the plasma. Comparing this result with the criterion obtained from 
the screening effect we conclude that, in an infinite plasma, multiple scattering is a more 
efficient mechanism for the dissociation of heavy mesons. 

The paper is organized as follows. In Sec. II, by analyzing the Dirac equation, we 
conclude that those photons with < I/clb and \kz\ ^1/0.3 are essential for the binding 
in the partonic language and get a Schrodinger-like equation in light-cone gauge in light- 
cone coordinates. We also discuss the typical time scale tb = i/Eb in a heavy meson. 
In Sec. Ill, we calculate the effective field induced by a fast moving charge in light-cone 
gauge up to the first order in ~ and use the uncertainty principle to estimate the 
criterion for the dissociation of fast moving heavy mesons. In Sec. IV, we give a parametric 
estimate for the dissociation of heavy mesons due to multiple scattering in an infinite hot 
quark-gluon plasma. In the Appendix, we compare the classical field calculated from 
classical electrodynamics and from QFT to illustrate the connection between the classical 
field and virtual photons. 
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II. THE DIRAC EQUATION FOR A FAST MOVING BOUND STATE 

Even though the correct mathematical treatment of a relativistic two-body system can be 
done by means of the Bethe-Salpeter equation, we choose to use the Dirac equation instead 
in order to get a simpler physical picture. This approximation is valid when one of the 
particles involved is much heavier than the other, in which case the field produced by the 
heavier particle is no longer a dynamical variable and can be treated as an external field. In 
this framework, we can investigate the general properties of the wave function of the lighter 
quark and determine under which conditions it will be bound. In order to be able to treat 
the system perturbatively we still have to assume both masses are much larger than Aqcd 
with one of the masses much greater than the other. Although the formal results are only 
valid in this case, we believe the parametric result should be the same for the case of equal 
masses. We assume that the plasma is in the deconfinement phase and in the perturbative 
regime. In this circumstance, we can treat the problem in the quark-gluon plasma in a 
similar way as that in a QED plasma. 

In the following we first investigate the role played by photons with momenta in different 
regions in the binding of a fast moving bound state and address the question of what is the 
approximate equation, appropriate for such a system, analogous to the Schrodinger equation 
for a bound state at rest. Then, we answer the question under which circumstances our 
analysis of the QED bound states applies to heavy mesons. We also give a brief illustration 
about the typical time scale tb = 'y/Es based on the perturbative definition of the wave 
function at the end of this section. 

Since the detailed screening effect for a moving bound state is too complicated to be 
solved analytically js], we try to simplify the problem and make it suitable for an intuitive 
understanding by analyzing the Dirac equation for a fast moving bound state and deter- 
mining which photons are essential for the binding. This will allow us to make appropriate 
approximations to the calculation of the effective field induced by the heavier particle in the 
presence of the plasma. With the effective field, we still need a Schrodinger-like equation 
in light-cone coordinates which manifests the uncertainty principle and simplifies the spin 
structure in order to enable us to use a physical analysis of fast moving bound states similar 



to that for bound states at rest in Ref. 



2j. In the following analysis, we assume that the 



lighter particle has a mass rria and a charge Ca while the heavier particle has a mass rrib 
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and a charge Cb- The whole analysis in the following applies to heavy mesons provided the 
quarks are heavy enough. 



A. Which photons are responsible for the binding? 

In the vacuum case the field A^^ is much simpler, in certain gauges, in the rest frame of 
the bound state. Therefore, Let us start with this simpler case and assume that we have 
already solved \l/on(p), the wave function in the rest frame of the bound state, which peaks at 
p = with a width Ap = l/as, and = En- By boosting to the lab frame, we obtain \l/n(p) 
which peaks at po± = with Ap± = I/qb, Poz = with Ap^ = 'j/aB and Pq = 7-En- By 

writing 

w^^f^;"!). (2) 

\X{P) J 

and inserting it into the Dirac equation in momentum space 

- ^4p) =eaj^, [4mn{p - k)] , (3) 
we get, in the chiral represent at ion [gI, 



P ■ (^X{p) - ma<^{p) \ f d^k A{k) ■ ax{p - k) 



p ■ a(p{p) — niaxip) j (2^)"^ I A{k) ■ a(p{p — k) 



(4) 



Given the above general properties of \l/„, we easily see that A^^{k) with < l/a^ and 
\kz\ ^i/clb gives the predominant contribution to the k integration in the right hand side of 
(j4]). As reviewed in the Appendix, in QED the classical field A'^{vkz, k) arises from virtual 
photons with momenta k in the wave-function of the heavier particle in the same gauge. 
Therefore, we conclude in the partonic language that those photons with < l/a^ and 
\kz\ ^l/dB are essential for the binding of a bound state moving at velocity v. 

This analysis allows us to give a qualitative picture of how the binding is affected by the 
presence of the plasma and under which circumstances the existence of bound states is not 
allowed. As will be seen in the next section, the presence of the plasma becomes manifest 
through an effective photon mass fi, which causes the corresponding screening effect. This 
effective mass provides a cut-off on the lower limit of the integration on the right hand side 
of (jll) and, therefore, determines if the photons responsible for the binding are still available. 
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Following this argument, the criterion for the dissociation of bound states at rest in a plasma 
is as — 

Back to the limit f ~ 1, the arguments stated above show that we can focus only on the 
contribution to the field from photons with ~ 'yk±. Before doing this, we need a simplified 
Dirac equation in light-cone coordinates which will enable us to use an uncertainty principle 
analysis for a fast moving bound state. 

B. A Schrodinger-like equation in light-cone gauge in light-cone coordinates 

In order to determine the existence of fast moving bound states in the presence of the 
plasma we would have to solve the Dirac equation corresponding to that system. For the 
present case, the Dirac equation is too complicated to be solved exactly. Nevertheless, we 
can get an approximate Schrodinger-like equation which will put us on familiar grounds to 
make order of magnitude estimates on the conditions under which bound states are allowed 
to exist. The major difficulty comes from the fact that in the rest frame of the charge the 
effective field is highly anisotropic in the limit f ~ 1[5]. In contrast, as showed in Sec. IIII CI 
in the reference frame with the plasma at rest, the anisotropy of the effective field in light- 
cone gauge is predominately due to the Lorentz contraction and can be easily handled when 
expressed in light-cone coordinates. This allows us to neglect the anisotropy of other kinds 
and greatly simplify the problem. By using light-cone coordinates and keeping terms up to 
(9(q;^) we are able to derive an equation in which, by means of the uncertainty principle, we 
can establish necessary conditions for the existence of a bound state. 

In coordinate space, the Dirac equation (jl]) takes the form 





from which we can get two second order differential equations 



{P 



CaA) ■ a{p - CaA) ■ a(p = mlip 
CaA) ■ a{p - CaA) ■ ax = mix 



(6) 



After some algebra we get 




{P 



CaAf + ea{B + iE) - a ip = mlip 
CaAf + ea{B -iE)-a x = mix 



(7) 
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In the following we only keep terms with expectation value up to (9(q;^) and neglect terms 
which are higher order in a. Under this assumption the second term can be ignored in Eq. 
(I7j) in both of the equations given there. We show a detailed calculation of the potential for 
the vacuum case in the Appendix which supports this statement. In the plasma the electric 
and magnetic fields are screened and are even weaker than in the vacuum case. Since the 
two equations are the same in this approximation, let us focus on one of them 



(8) 



Even though the dominant part of the potential is in the transverse components, the main 
contribution in the equation above comes from the A~ component since it is enhanced by a 
factor. In this way we get a Schrodinger-like equation in light-cone gauge in light-cone 
coordinates 



Pj + ml 
p+ 



+ eaA- 



(9) 



where = p'^ ± p^. Here we have assumed p± is of order ^ ~ arua and then we have 
dropped terms which are higher order in a. In this paper, we will not solve exactly, 
but instead we will use the uncertainty principle to estimate the existence of bound state 
solutions in the limit f ~ 1. Inspired by ([9]), we start with 



{p-)^^^7Z^ + eaA-{{r)), 



{p+) 



(10) 



where = - Ap+^J^ ~ ^/{p+^) 



2-fma 



2(27m„)2 



and we have 



taken a/ — 27ma. In the case that A (r) only depends on r = ^/x'j_ + 'y'^{x~y, that 
is, the system has a generalized rotational symmetry, we may expect 

if) 



m„ 



(11) 



with p = (p_L,Aj9+/7) and x = {x±,'-fx~). Except for the different definitions of the 3- 
components of x and p, the physical meaning of (fTTl) is exactly the same as that used in the 
uncertainty principle analysis from the Schrodinger equation. In the vacuum case, we have 



A- = l^^^i^ in light -cone gauge as calculated in the Appendix, and 



l{p > 



2m„ 2 



a 



(12) 
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2 

in the limit 7^1. By using the uncertainty principle we have 7 ~ ^ — " 



and p ~ i ~ ania- This is indeed consistent with the results obtained from boosted 
wave functions by keeping terms up to (9(q;^) in binding energy. Therefore, as in the limit 
t> ~ 0, we can get a Schrodinger-like equation in the limit f ~ 1 which allows us to use the 
uncertainty principle to estimate the properties of fast moving bound states. 

Before calculating the plasma effect on the effective potential, we give a quantitative 
estimate about how well the color Coulomb (perturbative) potential applies to the heavy 
mesons in the vacuum. Using the uncertainty principle with the Cornell confining potential 
V{r) = Kr - ^^], we have 



which tells us that if 



Knia 

P = — 5 h aCpma, (13) 

we may neglect the non-perturbative linear potential responsible for the confinement in V{r). 
In this case the discussion about the electromagnetic bound state given in this section is also 
valid for the heavy meson if we replace a with aCp- For the charm quark, ^^^^^g GeV^^ ~ 
l.SGeV^ with ml = ~ l.GGeV^ and a{mc) = 0.4. For the bottom quark, 

(^g^Ge1/2 ^ io.5Ge\/2 ^^^-^ ml = A.7^GeV^ - 22GeV^ and a{mb) = 0.2. Even though 
for the charm and the bottom, equation ( |T^ is not perfectly satisfied, as a parametric 
estimate in the following sections we can still take the binding energy Eb — a'^Cpma- 



C. The intrinsic time-scale in a bound state 

From a perturbative point of view, the interaction between the two quarks in a heavy 
meson occurs via interchange of gluons. It is possible to define a typical time within which 
we can neglect the interaction between them. This time scale will play an important role in 
establishing the appropriate criterion for dissociation due to multiple scattering. 

Let us start with the perturbative definition of the wave function^], 

^ r (15) 

- te^ j d% [SUxa - x'J4K)U 
where da{x'^) is the volume element of the closed 3-dimensional surface of a region of space 
time containing Xa, n'^{x'^) is the inward drawn unit normal vector of this surface at x'^, A'^ 
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^(t)>= * \'¥(to)> +S ^Y(t')> 



t' 



Y 

Eb 



T(to)> 



FIG. 1: The perturbative definition of a heavy meson. Black dots here represent kicks to quark a 
from the color field generated by quark b. lit — to < 'y/Es, the wave packet propagated by the free 
propagator represents the predominant contribution to the wave function at t. Approximately, the 
picture of the interaction between the two quarks in a heavy meson is that quark a is kicked once 
by the color field of quark b within a period At < tb = 'j/Eb to pull it back into the heavy meson. 

is calculated perturbatively by the one-gluon exchange approximation and we neglect the 
linear potential responsible for the confinement since the masses are assumed much greater 
than Aqcd- Choosing the integration surface on the first term on the right hand side of f|T5|) 
as the whole space at two different fixed times, we can relate the wave function at time ta 
to the values of the wave function at a previous time to- If '^t = ta — to <^ tb — 7^7^ — , 
the first term on the right-hand side of ( JTSl) gives the predominant contribution to the wave 
function at ta- Therefore, by conservation of probability, one can neglect the second term, 
that is, the interaction between these two quarks can be neglected. This is easy to see in 
light-cone coordinates, where the free propagator is 9| 

Sf{x) = I [u{p)u{p)e-'P-''e{x+) - v{p)v{p)e'''-^e{-x+)] , (16) 

with p~ = Ejjt^^ p = and p^ is integrated over the region < p+ < 00. If 

typically ^Ax"*" ~ " ^m™" — ^'^/'^b ^ 1, one can neglect it in the exponentials in flTB]) 
and, therefore, the integration of in the first term on the right-hand side of f|T5|) simply 
gives a (5-function to reproduce the wave function at ta in the limit 7^1. Otherwise, if 
^Ax~^ ~ " ^m'"" ~ ^^/'^B ^ 1, the first term on the right-hand side of f|T5l) is highly 
suppressed due to this big term in the exponentials in the integration and, therefore, 
the second term on the right-hand side of f[T^ contributes predominately by conservation 
of probability. Approximately, the picture of interaction between the two quarks in such a 
heavy meson is that the lighter quark interacts with the color field generated by the heavier 
quark once every interval At < tb = '^/Eb- We illustrate this in Fig. [TJ 
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III. DYNAMIC DEBYE SCREENING EFFECT FOR FAST MOVING BOUND 

STATES 

The next step following our analysis of the Dirac equation would be to calculate the 
appropriate potential for a fast moving particle to plug into equation ( ITTi) . However, this 
whole analysis ignores completely the effect of multiple scattering with the particles in the 
medium, which can modify strongly the wave function of the system. Although in general 
these two mechanisms interfere, we will treat them separately in order to determine which 
one is dominant and which one should be used to establish a criterion for existence of bound 
states in the plasma. This approximation relies on the fact that screening is a coherent 
effect involving a correlated motion of particles in the medium, unlike multiple scattering 
which is an incoherent effect caused by random kicks from uncorrelated scattering centers. 
In this section, we will give a detailed calculation for the hot QED plasma, and generalize 
the results to the quark-gluon plasma. 



A. The photon polarization vector within the HTL approximation 

First, let us calculate the screening effect on the field induced by the heavy particle, due 
to the presence of the plasma. This is done by calculating the retarded photon propagator 
in thermal field theory for different regions of momenta. Following the analysis presented 
in the previous section, we are mainly interested in those photons with momenta ~ 7fc±, 
which are essential for the binding of a fast moving bound state. In the following calculation 
of the retarded photon propagator, we keep terms only up to first order in kj_/kz ~ ^■ 
Moreover, we neglect the modification of A'^ due to the appearance of the lighter particle 
(for a detailed discussion, from the kinetic theory point of view, about the effective field 

in the rest frame of the particles and the influence of the appearance of another heavy 
particle see Ref. jsj). Even though the photon polarization vector within the hard thermal 
loop (HTL) approximation is well-known (say, 0, 111), in the following we still present 
some details of the calculation, which enables us to see how well our approximation is in the 
regions of momenta beyond the HTL approximation. 
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Using the notation in Ref. [10|, let us calculate the photon polarization vector [10|, lU] 



d^q Tr [riYii - m 



(17) 



with P^" = -g"^ + ^ - Pf, Pt^ = Sij - kkj^ Pt^ = pf = and cu„ = (2n + l)7rT. From 



Eq. ([T71) . we have 



and 



ttOO l2 

poo ^2^^ ' 



G = -l(F + n) 



(18) 
(19) 



with n = g^jjl^^ . After taking the trace and performing the frequency sums we get, 

1 1 



+ 1 



El-k-q 



E1E2 



n{E{) - fiiE^)) 



kO-Ei- E2 k^ + Ei + E2 
1 



k^ + Ei- E2 k^-Ei + E2 



(20) 



and 



El-k-q\ / 1 



1 



1 



+ 1 



E1E2 



{n{E,)-n{E2)) 



kO-E^- E2 k^ + Ei + E2 
1 



ko + E^- E2 k^-E^ + E2 



(21) 



where n{E) = -^^w^, E^ = \q\, and E2 = |g - k\. 

In the HTL approximation, that is, kz,k± <C T, we get the well-known results 



(Pq dn{q) 
(27r)3 dq 

A;0 



2 - — In 



2 - 



+ 1 



fco + /c cos 6 k^ — k cos ^ 



(22) 



k'" 



and 



n ~ 46^ 
^ -4e 



1 (1 - 2n(g)) 



(27r)3g 
2 / c?^? 1... 



(23) 



(27r)3g 



2n{q) = -2/i' 
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where we have dropped the T-independent divergent term and we have taken /x^ = 
When making use of the classical current induced by a fast moving particle, we introduce a 
delta function of the form S{k^ — vk^), so we can safely replace = vk^ in the calculation 
of the propagator. In our case kz — "ykx, by keeping only the first order of we have 

2-^ln^ U2A-^ln7^0, (24) 



P\ k f - 1 



and, therefore. 



G ~ jj'. (25) 



B. The photon polarization vector beyond the HTL approximation 

As explained in previous sections, photons with kz — 'y/aB play an important role in 
the calculation of bound state wave functions, and since we are interested in fast moving 
particles, kz might be comparable with or even larger than T in the limit f ^ 1. On the 
other hand, the transverse components of the photon momenta are not necessarily large. 
The relevant region for the transverse momenta is ~ I /as and in particular we would like 
to focus on the region l/a^ ~ a^^^T, which means we still have the condition k± -C T. In 
the following, we give a detailed calculation of the photon polarization vector in the limit 
k±<^T and kz > T. 

Let us calculate F and G for this region of momenta. Since we are considering the case 
I A; I ^ T, the predominant contributions to the integrals in fl^Ul) and fl^Tl) come from the 
separate regions with |g| ~Tor \q — k\ ~T. Both regions give the same contribution since 
the integrands are unchanged under q ^ k — q. Taking into account the contribution from 
both regions we get 

F = — 4e^ / ,^ ^ = —n' / dcosO- —. (26) 

p J (27r)3 q l-cos2^ p Jo 1 - cos2 ^ 

The angular integration has a coUinear divergence that is cutoff by the k± term neglected 
in the approximation above and by the mass of the constituents of the plasma (neglected in 
our calculation for high T). Therefore, as in the case kz, k± <^ T in (J2^ . this logarithmic 
divergence must be proportional to In 7, that is, 

F~ 7-2 In 7. (27) 
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Similarly, the predominant contributions of g-integration from the (separate) regions with 
|g| ~ T or \q — k\ ^ T gives us 

n = -2/x2. (28) 

Accordingly, 

Therefore, in the case kj_ <^ T, we have F ~ and G ~ /i^ in the limit ^ T as well as 
in the limit k^ <^ T. This justisfies F ^ and G ~ yU^ in the case fc^ <^ T in the following 
calculation even though the region with A;^ ~ T is difficult to evaluate analytically. 



C. The effective field in light-cone gauge 



To use (ITT]) for a fast moving bound state in the plasma, we need to evaluate the retarded 
photon propagator in light-cone gauge. By only keeping terms up to first order in k±/kz — 
we have 



(30) 



and 



1 kf^k'^ 



(31) 



where -D^^, defined in ( I5^ . is the vacuum light-cone gauge propagator. Using (1301) and 
f PT]) . we get the full photon propagator 



r]>'k'' + r]''k>' 



k'^ — jj? 



9' 



7] ■ k 



(32) 



The corresponding effective field is given by Maxwell's equations 

V — 1 k\ 1 ~ V 



A^'{k) = -iD^^^{k)j^ ^ 27Teb6{uj - vk,){ 



k"^ — fi"^' kz{k'^ — fJ^'^) ' k"^ — fi"^' 



and in coordinate space we have 



-fir 



A-(x) = ^2(l-.)7^, 



(33) 



(34) 
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and 



^.|^[e(a:-)-0(-a:-)]v. J dk.k,^^0 (35) 

- -r"-T/^a;_Li^i(/ia;_L) [6(x") - Q{-x~)] , 
47r x]_ '- 

where r = y'a;5_ + 72(z — ft)^ ~ \/ x']_ + 'y'^{x~)'^ and in the calculation of A±{x) we have 
only picked up the poles at = ±ie. 

In kinetic theory, the effective field A'^{k) is found to be strongly anisotropic for f ~ 1 in 
the rest frame of the charge [5], which seems to contradict the results in light-cone gauge. 
However, we shall see that after some gauge transformation and going back to the rest frame 
of the plasma, the anisotropy is suppressed as inverse powers of 7 in contrast with that due 
to the Lorentz contraction. In the limit ~ 1, the effective field A'^(k) in the rest frame of 

n 

the charge calculated in Ref. [5|] is 

1 1/1 1 



A'^{k) ~ 27reb6{uj 
A^(k) ~ —27iebS{uj) tan 6' cos 
A^{k) ~ — 27refc(5(co') tan^ sin 
A^{k) ~ 27rebS{u)tan'^e 



C0S2^ VA?2 +/i2 p 

1 1 

fc2 + ^2 1^2 



1 1 ^ ^''^ 



A;2 + fi'^ A;2 
1 



,A;2 + ^2 

with k^ = {uj, \k\ cos sin 6', \k\ sin sin 6*, \k\ cos 6'). Equ. fl36l) is the same as calculated in the 
above approximation F ~ and G ~ /i^ in covariant gauge. After the gauge transformation 
Af'{k) A^'{k) + k^'A{k) with 



we have 



A%k)-27reb5{ij) ^ _ . . 

LA;2 + ^2 cos2^ VA;2 + /i2 k^ , 

i±(A;)~0, (38) 



1 
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Back to the rest frame of the plasma, it gives 

A\k) ~ -2nebS{uj - vk,) ^ 

/e^ — jjL^ 

A^{k) ~ 0, (39) 

showing that the anisotropy is suppressed as inverse powers of 7 in the rest frame of the 
plasma. We have shown that the effective field of a fast moving charge has a smooth limit as 
f — *■ 1. Moreover, by using light-cone gauge in the rest frame of the plasma, one can obtain 
an effective field in which the anisotropy is predominantly due to Lorentz contraction. This 
allows us to use ffTTl) to obtain a criterion for the dissociation of fast moving heavy mesons 
due to dynamic Debye screening. 

D. Dissociation of heavy mesons due to dynamic Debye screening 

In the HTL approximation, one can obtain the gluon self-energy simply by taking /i^ = 
^g'^T'^{Nc + ^) in fl2^ and fl25|) 11|. where Nf is the number of massless flavors of quarks 



in the quark-gluon plasma. Since this result does not depend on the gauge-fixing [12], by 
replacing et with —g, we get the effective gluon field in light-cone gauge 

A-(x) = -f 2(l-t;)7^, (40) 
Att r 

and 

A^{x) ~ -^^f,x^K,{fix^) [Q{x-) - Qi-x-)] . (41) 

Now, we are ready to discuss the criterion for the dissociation of fast moving heavy mesons 
in the quark-gluon plasma by using the uncertainty principle in light-cone gauge. Inserting 
fl34D into (fllD, we have 



^{p)-o — + ^-2(1-^)7 — + ^ , (42) 

2ma 2 r 2ma 2 r 

in the limit f ~ 1. Assuming r = 1/p and minimizing fH21) . we have 

^ -x(l + x)e-^ (43) 



aCpma 

with X = jj/p. Eq. (H3!) has solutions only if /ia^ = < 0.84 (see [^). Therefore, we 

can use the uncertainty principle argument for the dissociation of fast moving heavy mesons 
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in the plasma in the same way as it was used for bound states at rest The corresponding 
criterion for the dissociation of fast moving heavy mesons, based on a screening analysis, is 

ttB ^ (44) 

Note, the anisotropy in the effective field seen in the rest frame of the charge could imply a 
more efficient screening than that in the case with the charge almost at rest with the plasma, 
but the modification of this criterion on the right hand of (jH]) should not depend on inverse 
powers of 7. 



IV. DISSOCIATION OF HEAVY MESONS DUE TO MULTIPLE SCATTERING 



In this section, we give a parametric estimate of the criterion for the dissociation of 
heavy mesons due to multiple scattering in terms of the saturation momentum Qs, which is 
a characteristic property of any QCD mediaQ, [l^. Quantitatively, the physical meaning 
of Qs is that the gluon distribution of the target is dense as seen by probes with virtuality 
q± Qs, but it is dilute as seen by probes with high virtuality. When a heavy meson with 
a size travels in the quark-gluon plasma, l/as naturally plays the role of the virtuality 
q±. If ttB ^ the meson will break up, that is, the medium looks opaque to the meson. 

This picture is confirmed by detailed calculations in Ref. 1^. 

This picture can also be justified by the argument of the transverse momentum broadening 
of the two quarks in a heavy meson. As showed in Fig. [1], within a period tb, the two particles 
in a heavy meson propagate in the medium like two free quarks and pick up the transverse 
momentum broadening (Ap^), which is equal to the saturation momentum squared Q^jl^. 
If (Ap^) = ^ l/'^B) that is, the two quarks pick up transverse momenta greater than the 
typical momentum in a bound state, the meson will break up. Therefore, the criterion for 
the dissociation of a heavy meson due to multiple scattering in a medium is parametrically 
aB — 1/Qs- This picture applies to cold matter as well as to hot matter |l^. In the following, 
we shall give a parametric analysis of the criterion for the dissociation of heavy mesons in a 
hot quark-gluon plasma in the case that the successive scatterings are uncorrelated. 



In a hot quark-gluon plasma modeled 
quark, Q^ has the following simple form 



)y uncorrelated scattering centers, for a fast moving 



16| 



Q'i ~ In — = Lpa/i^ In — ~ o'CpiN, + -^)T^L In (45) 
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where we have taken the number density p ~ [Nc + ^)T^, a = ^^C^ and /i^ = ^{Nc + 
^)g'^T'^ ~ a{Nc + ^)T^ = a^ffT"^- In a finite plasma with a length L < Tb, the criterion 
for the dissociation of heavy mesons is a ^Cf{N^ + ^) In ^T^L ^ ^. In an infinite quark- 
gluon plasma, the time scale tb = -3- — -^jr- plays the role of the length L in the definition of 



Qs, and our criterion for dissociation of heavy mesons — 1/a^ gives ~ ^ 



for the quark-gluon plasma, in contrast with that due to the screening effect as — - — —r — ■ 

"?//^ 

If the successive scatterings between the mesons and the plasma constituents are essen- 
tially independent of each other, the criterion as — ^ — is also true for the disso- 

elation of heavy mesons almost at rest with the plasma. Quantitatively, in a hot plasma this 
means the collision time Tc — \ = — — } rp ^ 1/yU, which is equivalent to 1 ~ %a ^ ^ «• 
In this case, non-relativistic quarks pick momentum broadening symmetrically in each direc- 
tion due to the uncorrelated random kicks from plasma constituents in the same way as the 
transverse momentum broadening of relativistic quarks. Therefore, we would expect that 
the momentum broadening of non-relativistic quarks should have a similar form as (H5|) . 



Since uncorre 



ated multiple scattering implies a <^ 1, we can use the leading-log 



approximation 171] in our discussion about the criterion for the dissociation of non-relativistic 



heavy mesons. The mean-squared momentum transfer per unit time between a non 



relativistic heavy quark and the plasma is calculated in Ref. 
the leading log terms, we have 



and by keeping only 



I (lApf) ^ ^a'C,{N, + ^) In (46) 
Since Eq. P6|) is time-independent. 



{{Apr) ^ ^a'CAN. + ^) In ^T^At. (47) 

Taking At ~ ~ and ((Ap)^) ~ I/c^b, we get ~ j—. Therefore, we 

conclude that as — — is a parametric criterion for the dissociation of a heavy 

meson due to uncorrelated multiple scattering in an infinite plasma. This criterion was also 
obtained in Ref. ^19] in the case of heavy quarkonia at rest with the plasm a by effective 
field theory techniques, which was first obtained (without the logarithm) in 20|. At high 
energies, heavy mesons can dissociate even at much lower temperature T before dynamic 
Debye screening effect plays an important role. 
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APPENDIX: The Classical Field A'^^ in Coulomb and Light-cone Gauges 

In this appendix, we compare the field A^^ for a charge Cb moving with a velocity v 
along the z direction in classical electrodynamics and the classical field A'^i for a charged 
particle with p'^ = p^{l,0±,v) and charge Cb in the partonic picture in QED. In classical 
electrodynamics, the field is calculated using Maxwell's equations with a classical current. 
In momentum space, we have 

A>'{k) = -iD^'{k)Uk), (48) 

and 

3^{k) = 2'Keb5{u - vkz) (1, Ox, v) . (49) 

(i) Coulomb gauge v ' ^ = 

The photon propagator in Coulomb gauge is 



^ 



with = 6ij — kikj. The classical electromagnetic field is 



yi3 



A^{k) = -tD^J{k)Uk) = 27Teb6{uj - vkz){^, -^)- (51) 

k^ fc 



In the limit k, ^ k\ and f ~ 1, we have 



rv^ rv 

and 



A'^ik) ^ 2neb5iu - vkz){y^, 1), (53) 



which has the same transverse components as (l56l) . 
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|P)= 



P P-k 

tk,A 



FIG. 2: The dressed wave function of a charged particle 



(ii) Light-cone gauge A'^ = 
The photon propagator is 



Tj ■ k 



and 



A^{k) = -iD>l'cik)Uk) = 2net6iuj - vk,){- 



1 — V k\ 1 — f , 



where rj'^ = 0, 0, —1). In the limit f ~ 1, we have 



A^'ik) 2nebS{uj - vk,){0, 



0), 



(54) 



(55) 



(56) 



kzk\ 

which is the same as the Weizsacher- Williams field calculated in the light-cone wave function 
of particle M„ the parto„,c picture Q. 

In QED, the field is quantized and virtual photons appear in the wave function of a 
charged particle to give rise to the classical field. In this appendix, we use the notations in 
Ref. ||6|] and the classical field is defined byfl 



{p'\A^'{x) \p) 



(57) 



2Ep.(27r)3 

where A^{x) is the quantized photon field and \p) is the dressed wave function of the charge 
particle as showed in Fig. [2] 

d^k 

\=± ' 
with 

{p' -k;k,X\Hj\p) 

(59) 



x{k) \p-k;k,X) 



(58) 



2k\27r)'2E._,;^x{k)6{p-P 
= eb{2Txf6{p-p 



p^-{p- kf - fco 
u{p - k)^\{k)u{j)) 



pO-{p- kf - ko ■ 
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Inserting (158|) into (157|) . we have 

^cii-) = E / (03^^'-"^a(^)^a(A:) = / (60) 

Here, we only keep the positive-energy part of A^^. 

The classical field A'^i in light-cone gauge is well-known jl^ . And we shall calculate it in 
Coulomb gauge as another example to illustrate the correspondence between A^^ and A'^i. In 
Coulomb gauge, 

E <^A* = Pf- (61) 

Assuming , we have 

^.{k) ^ -'-^S^ ^ -e,J^^ = e.^i^, (62) 
^ ' 2 k^{l-vk,) k\l-v^kl) k^ ' ^ ^ 



and 



AMA:) - e,-^^^ ^ e,^, (63) 



with k^ = v'^k'l — k"^, which, if multiplied by 2tt6{uj — vk^), is the same as the vector potential 
in ( |5T1) . Therefore, both calculations in Coulomb and light-cone gauges illustrate the fact 
that the field A^{uj,k) with uj = vk^ in a certain gauge in classical electrodynamics arises 
from the virtual photon with momenta k'^ = {\k\,k) in the same gauge in QED. 

Since the Dirac equation is more convenient to solve in coordinate space, we shall evaluate 
the field y4'^(x) by the Fourier transformation of flSTl) and fl55l) . In Coulomb gauge, A^^lx) is 

d^k 



/ .2 X X (64) 



1 / 7 

6^ I V A(x) V I ~h A(x 

Attt ' dx\dz ' \ Attv dz^ 
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where r = ^Jx\ + {z — vty, r = ^/ + 7^(z — vt^ and 



d^k e ^^'^ f cPk e ^^'^ 



(2vr)3 + ^-2^2) y (2vr)3 ' [A:i + (r/ + (1 - vh-')klY 

^ r dk±dk,k± e'''^^'-''*\Joik±x±) 

"-V I 77777^5 77T~7Z~77~Vu^72u:2T2 



'J (27r)2 [A;i + (r^ + (l-r/)7-2)A:2] 



^ rf?7 + (1 — ?7)7~^)x^ + (-2 — vty 
Stt + (1 — ?7)7~^ 

]^ r — 7"^f — (z — ft)ArcTanh (^z^) + — ft)ArcTanh ( 



7 -^r 
V z—vt 



471- _ 1 

Inserting (165|) into (16^ . we obtain 



^.(,) = £i (i, _ 1 _ i), i (i _ i 1 ) ) . (66) 



In light-cone gauge, in coordinate space, 



(27r)3 [kl + 7-2fc2] A;, 
/" d'^k _iu.ri. 1 



V- 



(27r)3 [A;i + 7-2^2] k, 
dkjk±ki_ jkAz-vt) Mk±x±) 



(27r)2 [kl + -f-^k^] K 

/dk Qikz{z~vt) I 
(2vr)2 k. M^k^j-^xl) 

-^VxArcSinh^ii^ 

Att x± 

Cb x± -f{z - vt) 

4:71 X^l f ' 



and 



In the limit v 



^({l-vhiV-^. 

47r \ r Xj_ r r 



A^{x) = ^ ( (1 - v),-^ - ii±^, _(i _ ) . (68) 



(67) 



A^{x) = ^ [e{x-) - e{-x~)] , 0^ , (69) 

which is the same as A'^i{x) in Ref. 14 1. 



^ Note, in Sec. II and Sec. Ill, this is defined as r. 
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It is interesting to notice the big difference of the role played by different components 
of the classical field A'^{x) in light-cone gauge in scattering processes at high energies and 
in fast moving bound systems. In a scattering process at high energies, if we take the 
right-mover as a classical current, since the other particle involved is a left-mover, only the 
transverse cornponents A± of the right-mover give a dominant contribution to the amplitude 
of the process IJ]. On the other hand, in a fast moving bound state, since the two particles 
involved are both, say, right-movers, it is A~ that plays a more important role since it apears 
in the Dirac equation via the product with p+ to contribute a term of order while the 
contribution from A± with the product of p± is of order a^. 
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